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I. INTRODUCTION 

One of the primary tasks of materials scientists 
is the preparation of solids with controlled shapes, 
compositions and structures. These characteristics in 
turn determine the properties of the solid material. The 
degree of our ability to exercise control over them 
ultimately limits the performance of the material for 
its intended application, whether it be structural, electronic, 
magnetic, optical, chemical or biological. Consequently 
studies of materials processing are highly important areas 
of research which- involve interdisciplinary combinations 
of, such diverse fields as solid state physics, surface 
physics and chemistry, fluid dynamics, heat and mass transfer, 
and high temperature thermodynamics and physical chemistry. 
Although there are many important processes such as metals 
casting and annealing, ceramic sintering, glass melting 
and homogenization, chemical vapor deposition, crystal 
growth, oxidation, and diffusion, this paper will consider 
only crystal growth from melts and solidification processes 
since these are uniquely affected by a reduced gravity. The 
basic configurations and problems associated with crystal 
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growth from melts will be outlined on earth as well as 
in space. The role of free liquid surfaces in crystal 
growth processes will be discussed for various liquid 
shapes and the relevance of the current state of the 
science of liquid drops and bubbles in liquids to these 
processes will be demonstrated. 

IX, CRYSTAL GROWTH AND SOLIDIFICATION PROCESSES ON EARTH 
Melt crystal growth techniques can be classified 
in terms of the degree of confinement of the liquid by 
an external container as shown in Fig. 1. In the full 
confinement methods, the growth interface contacts the 
container over some region as shown for the horizontal 
and vertical Bridgman techniques. ^ The main problems 
associated with this contact are direct container contamina- 
tion of the melt close to the interface itself and the 
thermal stresses occurring between the container and the 
grown crystal. In the partial confinement methods, the 
growth interface does not touch the melt container as shown 
for the Czochralski, flux and hydrothermal growth techniques. ^ 
Although contamination from the container does frequently occur 
in these methods, it can often be reduced through the use 
of special liners inside the structural support of a rigid 
container. The Czochralski technique is the most widely 
used method for growth of congruently melting materials 
by direct solidification. The growth interface is constrained 
by appropriate temperature gradients and the meniscus of the 
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free melt surface. Included here are the shape defining 
techniques such as the Stepanov method and edge-defined/ 
film-fed method, both of which will be described later. 

The various flux techniques are basically isothermal and 
are used for materials that are incongruently melting or 
decompose at higher temperatures. The interface position 
and growth rate are constrained by the compositional 
gradients in the melt and the degree of supersaturation of 
the melt respectively. For some fluxes, the free surface 
results in undesirable vaporization of some constituents. 
Hydrothermal growth produces crystals from a supersaturated 
fluid near or above its critical temperature. The growth 
rate is controlled by a concentration gradient induced by 
a temperature difference between the nutrient and growth 
regions. 

Crystal growth methods from unconfined or container- 
less melts involve both drops and zones as shown in Fig. 1. 
Constant volume sessile drops which are fed externally 
with "raw material” are used for both flux growth by the 
vapor-liquid-solid (VLS) technique and congruent solidifica- 
tion by the Verneuil technique. Similarly, liquid zones, 
where the melt surface is constrained by surface tension 
assisted in some cases by electromagnetic levitation, have 
been used for both congruent melting (floating zone, silver 
boat methods) and fluxed melt systems (temperature gradient 
zone melting or travelling solvent methods) . None of these 
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techniques has been widely used on earth because they are 
complicated to control and because they require rather 
high temperature gradient s due to t h e^I^ara ts impo s ed on 
the liquid volume size by hydrostatic pressure considera- 
tions, The high temperature gradients cause extensive 
dislocation generation in the grown crystals and also give 
rise to highly nonuniform and unstable thermal convection 
in the melt which in turn leads to undesirable compositional 
variations in the crystal through its influence on mass 
transport at the growth interface. The removal of gravita- 
tional constraints however leads us to reconsider the 
detailed behavior of unconfined liquids and the influence 
of this behavior on crystal growth processes in a later 
section. 

The use of bubbles in materials processing on 
earth has been mainly as a means of providing large surface 
areas for gas-liquid reactions in extractive metallurgical 
processes and for stirring purposes in glassraaking. In 
reduced gravity environments these techniques cannot be 
used because it would be difficult to remove the bubbles. 

On the other hand several novel uses of bubbles in materials 
processing are possible in reduced gravity which will be 
discussed in the next section. 

III. CRYSTAL GROWTH AND SOLIDIFICATION PROCESSES IN SPACE 

Some possible configurations for containerless 
crystal growth in space are shown in Fig. 2. A division 
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can be made Into single, double and multiple solid access 
into the liquid drop. 

In the single access techniques, the drop is 

solidified onto a crystal seed so that as the seed is 

withdrawn , the drop position must be maintained. There 

are two main techniques for holding the drop position; by 

12) 

external fields such as electromagnetic v ' and stationary 
(3 4 ) 

acoustic ' ' fields, and by the use of dies such as used 
in the Stepanov^ and edge-defined, film-fed techniques. 
The latter techniques may be used to grow noncircular 
cross-sectional shapes but reintroduce contamination 
problems. The use of external positioning fields possesses 
several complications of which the most important are 
uncontrolled convection, drop oscillations and what heating 
method is to be used. In the case of acoustic field 
positioning, drop oscillations (and associated fluid flow) 
may be present, optical image heating methods must be used 
to avoid distortion of the sound fields and the energy 
impedance mismatch which gives rise to the positioning force 
becomes weaker in high temperature gases. In the case of 
electromagnetic field positioning, only liquids with a 
reasonable electrical conductivity (to 1 mho cm”^) may be 
used and uncontrolled electromagnetic stirring may be present 
at higher power levels if the field is used to simultaneously 
heat as well as position the drop. Both positioning field 

arrangements will soon be tested in the space environment 

(8 ) 

achievable in sounding rockets. 
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In all single access drop techniques, control 
of the meniscus shape at the freezing growth interface 
will be required for ease of growth-rate and crystal shape 
control. The role of the meniscus shape will be discussed 
in the next section but it is desirable that the liquid 
and solid surfaces meet at an angle of 180°. 

The need for positioning fields may be eliminated 

by the use of double access liquid drop techniques where 

a melting rod of the same material as the solidifying 

crystal is introduced into the drop. The simplest concept 

is the cylindrical floating zone technique where the 

end rods are of equal diameter, in this case, the maximum 

stable liquid zone length is equal to the circumference 

( 12 ) 

regardless of the surface tension. 7 Although this zone 
length is longer than that possible in earth gravity, a 
problem arises because the temperature distribution in the 
absence of thermal convection is not favorable. As shown 
in Fig. 3 (b) , the interface shapes will not be perpendicular 
to the crystal axis as shown in Fig. 2 but will rather 
follow the freezing point isotherm resulting from thermal 
conduction. The simplest remedy to this problem will be to 
use unequal and diameters as typified by the pedestal 
technique in Fig. 2. This geometry will allow heaters 
to be used which possess axial as well as radial temperature 
gradients. If zone lengths longer than the crystal circum- 
ference are desired, then floating drop shapes may be used 
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as shown in Fig. 2. However independent meniscus shape 
control is not possible for such drop shapes. 

Multiple solid access in the liquid drop 
allows for compositional control of the molten zone 
during growth. In such a system, complex solid-solution 
systems could be grown at constant solid composition 
by continuously changing the liquid composition. However 
such a technique represents a degree of difficulty not 
presently attempted on earth with the exception of the 
feed- in, pull-out and floating crucible techniques. ^ 

The role of bubbles in liquid drCps has been 

( 13 ) 

shown in recent Skylab experiments ' ' to effectively 

I 

"stiffen” the outer liquid surface so that vibrational 
oscillation amplitudes are considerably diminished. 

Consequently any of the techniques shown in Fig. 2 may 
benefit by using liquid Shells instead of liquid drops 
if drop oscillations are a problem. 

Another area where bubbles may serve a useful 
function is in the production of foam solids - a new class 
of materials with potentially useful strength to weight ratios. 
IV. LIQUID DROP STATICS AND DYNAMICS 

In this section, we consider four major phenomena 
relevant to the space processing of molten materials; the 
degree to which a contained liquid adopts its equilibrium 
shape in zero gravity, the stable shapes of axisyrometric 
menisci at crystal-melt boundaries, the rotational stability 
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of liquid drop and zone surfaces, and the vibration- 

oscillation characteristics of liquid drops and zones* 

These phenomena are all basic to our understanding of 

container less handling of melts for crystal growth purposes. 

(a) Behavior of Confined Liquids at Zero Gravity 

There are two important factors which determine 

the equilibrium shape of a liquid in a container under 

weightless conditions? the filling factor and the liquid- 

container contact angle. This problem has been solved 

(14) 

analytically for a spherical container by Zenkevich ' 
and his results are shown in Fig. 4. The transitions to 
the various morphological forms are given by the two 
curves A and B which denote the contact angle, 0, at which 
the liquid transforms from its earthbound position at 
an initial height, h, in a container of radius R. Below 
curve A, the liquid uniformly wets the container while 
above curve B, the liquid drop does not contact the container 
at all. A similar set of conditions applies to liquids 
contained in cylindrical tubes although no analytical 
solutions are available as yet. These results show that 
completely confined liquids may adopt a variety of shapes 
which may not be desirable from the standpoint of crystal 
growth and solidification, particularly in the vicinity 
of the transition points of Fig. 4. 

Cbl Static Meniscus Behavior at Zero Gravity 

The behavior of the meniscus shape during 
Czochralski growth on earth has been discussed thoroughly 
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by Geist and Grosse. The stabilities of such menisci 

C 16 V 

have been dealt with by Huh and Scriven as well as 
Padday and Pitts. Basically the situation in 
Czochralski growth is shorn in Fig. 5 where a solid 
cylinder is immersed in the free liquid-surface at a 
contact angle of zero. The preferred meniscus height, 
h 2 , is shown in Fig. 5(b) since the freezing interface 
will tend to remain at constant diameter for small non- 
uniform displacements of the crystal pulling rate, V Q , 
or freezing point isotherm. The stabilities of menisci at 
zero gravity for the free surface, single access drop 
methods of Fig. 2 have not been investigated as yet but 
are going to depend quite sensitively on the configurations 
of the force fields used to constrain the liquid drop 
position. Some anticipated meniscus shapes at zero gravity 
are shown in Fig. 6. The maximum height, h, of the crystal 
interface above the original drop surface is determined 
primarily by the radial departure, R, allowed at the bottom 
of the meniscus. If R is constrained to be close to the 
crystal radius, R c , then the height h cannot exceed 2 irR 
without the meniscus becoming unstable. Otherwise longer 
liquid meniscus columns are allowed according to the 
approximate relation: 

l \ n3T % 2 itIR 2 -CS-R c ) 2 :! 1/2 


( 1 ) 
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These considerations apply only if the liquid column 

base radius, R, is substantially smaller than the liquid 

drop radius. Otherwise the constant volume condition of 

the drop or zone will enter into the stability picture 

as pressure changes in addition to those caused by the 

liquid surface curvature. A treatment of the mathematical 

stability of liquid surfaces of revolution under such 

(18) 

conditions has been given by Gillette and Dyson for 
equal end radii. A summary of their calculated stability 
limits is shown in Fig. 7 for liquid zones of volume V, 
length 1, radius R and diameter D. The limits for unequal 
end diameters and differential volumes have not been 
calculated or observed as yet. 

(c) Dynamic Behavior of Liquid Surfaces at Zero Gravity 
The behavior of free liquid surfaces under 
rotational and vibrational conditions is now examined. 

(i) Rotational Stability 

In crystal growth and solidification processes, 
rotation is usually required in order to provide radially 
symmetrical temperature distributions since such symmetry 
may not be present in the heating source and perfect alignment 
of the crystal axis with the thermal center of symmetry is 
difficult to achieve. For purposes of compositional 
homogeneity however, no convection is desired in the melt. 
Consequently only uniform "solid-body" rotation of the liquid 
Cequal isorotation in the case of double access drop growth 
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techniques 1 may be considered and the absolute rotation velocity, 

- 2 

u, must be small enough so that the centrifugal acceleration, m r, 
does not cause density gradient driven flows (rotation 
of a 10 cm radius drop at 5 rpm results in centrifugal 
accelerations of 2.5 * 1Q~ 2 g) . Nevertheless such rotations 
will influence the surface stability criteria discussed 
in the previous section. 

The rotational stability criteria for liquid 

(195 

spheres was first considered by Rayleigh and more 
completely by Chandrasekhar.^" 20 * The deformed shape 
considered was that of depression at the poles and expansion 
at the equator as shown in Fig. 8 for various values of the 
dimensionless rotation rate, pQ^R 2 /a (where p is the density 
and a the surface tension of the liquid drop of initial 
radius, R, being rotated at Cl rad/ sec) . A water drop of 
radius 10 cm would become toroidal at a rotation rate of 
12 rpm. However recent crude experiments performed on 
Skylab with rotating liquid drops indicated that rather 
than the expected toroidal deformation shown in Fig. 9(a), 
a dumb-bell deformation (Fig. 9(b)) was produced when even 
small nonaxial perturbations are introduced into the rotating 
drop. The growth rate of such nonaxial perturbations is 
considerably reduced as the viscosity of the liquid drop 
increases. Such instabilities have not yet been analyzed 
theoretically . 
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The rotational stability criteria for cylindrical 
liquid columns has been analyzed theoretically by Hocking 
and Michael ^ 21 * , Hocking ^ 22 * , Gillis ^ 23 ^ , and Ross ^ 24 ^ . 

For axisymmetric disturbances, the maximum zone length is 
reduced below that for a static zone: 


L . 2 irR 


C2) 


This criterion has been tested experimentally in a Plateau 

( 12 } 

simulation system by Carruthers and Grasso v ' and in real 

(25) 

zero gravity on Skylab. These data are summarized in 

Fig. 10 for equal isorotations and compared to Eq. (2) . 

Discrepancies in the Plateau results occur because the 

actual zone rotation velocity is somewhat smaller than the 

end-member rotation velocity because of the viscous drug 

effects exerted by the viscous outer fluid containing the 
( 2 ) 

zone. v ' The Skylab data also shows a deviation from the 
theoretical curve - becoming unstable at lower rotation rates 
than expected. The reason for this behavior is that the 
failure mode in the Skylab experiments was nonaxisymmetric 
as shown in Fig. 11 lb) . it was possible to obtain axisymmetric 
failure modes of the type shown in Fig . 11 (a) by replacing 
the water zone in the Skylab experiments with a soap solution/ 
air foam which greatly increased the effective viscosity* 

The nonaxisymmetric failures clearly arise from disturbances 
created by misalignment of the rotation axes and nonparallelism 
of the discs at the ends of the zone. However viscosity 
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plays a very important and unsuspected role in determining 
which types of disturbances are amplified. A great deal 
of theoretical work is required here to analyze the 
nonaxisymmetric failure modes for both rotating liquid 
spheres and cylinder. 

tiil Vibrational Behavior of Free Liquid Surfaces 

In crystal growth processes, the time dependent 
deformation of free liquid surfaces is undesirable because 
the meniscus shape at the freezing interface will change so 
as to cause crystal diameter perturbations and also because 
internal fluid flows will develop which influence the growth 
segregation behavior. There are two separate liquid shapes 
which can be considered here as was done for rotational 
behavior; spherical drops and cylindrical zones. We are 
interested in the natural resonant oscillation frequencies, 
which depend primarily on surface tension, and the decay 
rates, which depend mostly on viscosity. Knowledge of 
this behavior will permit apparatus designers to avoid 
vibrations which cause resonant, large amplitude oscillations 
in the liquid volumes used in space processing. 

Lamb first calculated the oscillation frequencies 


of spherical liquid drops by assuming the viscosity was 
C261 


zero as: 


•qa u+1) g-ii (&+2)*| 1/2 

R 3 Cp o a+Pi U+l) ) J 


8 * 


(3) 
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where A is an integer greater than unity and p Q and 
are the densities of the fluids outside and inside the 
drop respectively. The role of viscosity was determined 
analytically by Miller and Scriven v ' who found that the 
time dependent frequency, 6, varied according to: 

6 * g R + i (4) 

where is the oscillation frequency and S R the decay 
rate and both are viscosity dependent: 

(2 i+li 2 (e*y i P 0 P i P 0 ) 1/2 
Sr 2/1 Rr[ (p i P i ) 1/2 + Cp 0 p 0 ) 1/2 J 

C2 i+l)I2U 2 -l)y. 2 p.+2 £(A+2)u 2 p +y u. (p. CA+2)-p (i-1))] 

. Jm A V W 4. X Q /P » 

2 R 2 ri(u i p i ) 1/2 + Cv 0 p 0 ) 1/2 ] 2 

e . B . . (2 (6) 
2/1 Rrr(u i p i ) i/2 +(u o p o } 1/2 j 

Here and u Q are the dynamic viscosities of the inner and 
outer fluids respectively and 

r - p Q Jl + p i fi+l) C7) 

Experimental verification of these relationships 

can be found in the Plateau simulation work of Valentine 
(28 ) 

et. ale where drops of a carbon tetrachloride and benzene 
solution with a density of 1 gm cm"^ were suspended in water 
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and oscillated by coalescsnse with small cyclohexanol 

( 29 ) 

drops. More recent experiments using sessile drops 
of the same organic solution suspended in water are shown 
in Fig. 12. The drops were deployed by simply filling 
a glass laboratory funnel so as to obtain the desired 
size (Fig. 12 Cal and Cel 1 and then mechanically vibrated 
(Fig. 12(b) and (d)). The data are plotted in Fig. 13 
and show that although the decay rates agree well with 
Eq. (5 ) , the oscillation frequencies appear to lie between 
3* and Sj. The contact of the carbon tetrachloride/benzene 
drop to the glass tube may be- responsible for this 
discrepancy although the contact angle is 130° so that 
the drop does not wet tube itself. 

More recent experimental data has been obtained 
with liquid drop experiments performed for the author on 
Skylab mission XV by scientist-astronaut E. Gibson. In 
these experiments sessile water drops were deployed on a 
circular disc which was then oscillated along an axis 
perpendicular to its face. The data are shown on Fig. 14 
together with the theoretical predictions of Eqs. (3) and 
(5) . One free drop experiment performed by astronaut 
W. Pogue on the same mission is included for comparison. 

In all cases f agreement with theory is excellent for both 
free and sessile drops. 

The vibration behavior of cylindrical liquid zones 
was also studied extensively on Skylab IV. Preliminary 
data for the oscillation frequency of longitudinal modes 
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(vibration along the zone axis) are shown in Fig. 15 
together with the theoretical prediction (solid curves) 


by Lamb 


(26) 


for an inviscid fluids 
1/2 

l o <rr 0 I 




2 tr 


L Lp 


( 8 ) 


where L is the zone length. The frequency increase at 
longer zone lengths is not presently understood but may 
be a result of the viscous flow effects in the zone. The 
theoretical vibration analysis for zones should be redone 
including these effects as was done for spherical liquid 
drops (Figs. 13 and 14). 

V. SUMMARY 

The unique environment of space allows the 
containerless handling of melts to be used for controlled 
crystal growth and solidification processes. Some of the 
techniques currently used on earth together with their 
problems have been reviewed. Possible configurations 
for containerless methods in a space environment have been 
discussed. Such techniques require a detailed understanding 
of the behavior of free liquid surfaces under both static and 
dynamic conditions. The current state of knowledge in the 
science of liquid drops and bubbles relevant to containerless 
processing has been discussed together with the results of 
recent Plateau simulation and Sky lab experiments. Many areas 
which require further theoretical and experimental work have 
been identified as a result of this preliminary work. 
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BRIDGMAN VERTICAL 

BRIDGMAN 

FULL CONFINEMENT 



CZOCHRALSKI FLUX HYDROTHERMAL 

PARTIAL CONFINEMENT 



' VLS LEVITATED FLOATING LEVITATED ZONE 

VERNEUIL DROP ZONE (SILVER BOAT) 


UNCONFINED MELT TECHNIQUES 


Fig. 1. Schematic classification of iaelt crystal growth processes on earth 
in terms of the degree of liquid confinement by a container. 





CZOCHRALSKI STEPANOV EDGE-DEFINED 

FILM-FED 

SINGLE ACCESS DROP TECHNIQUES 



FLOATING ZONE PEDESTAL 

DOUBLE ACCESS DROP TECHNIQUES 




Fig. 2. 

Schematic classification of 
melt crystal growth techniques 
in space in terms of solid 
access to the liquid drop. 


MULTIPLE ACCESS DROP TECHNIQUES 
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(b) 


interface profiles in zones for (a) 
convection in the melt and (b) zero 
conduction occurring in the melt. 


earth gravity 
gravity with 



CONTACT ANGLE 
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Fig* 5* Menicus configurations important at the crystalmelt interface during 
Czochralskl growth on earth; (a) necking instability* (b) desired 
configuration* (c) crystal diameter control difficult* 
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Fig, 7. Stability region of the surface of cylindrical liquid zones of 

various shapes contained between equal diameter aid discs - from 
Gillette and Dyson. (18) 



186 



Fig, 8. Calculated shape deformations of an initially spherical liquid drop 

of radius, R, rotating at increasingly fast rotation rates. Rotation 
axis is the z-axis and r is the equatorial radius - from Chandrasekhar. ( 
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Flg« 10* Comparison of the rotational stability of cylindrical liquid rones in 
Plateau s imu lation and Skylab experiments with theory for axisymmetric 
modes. 
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Fig* 11* Forms of rotating liquid zone surfaces: (a) axisymmetrlc mode observed 
in presence of highly viscous liquids, (b) nonaxlsymmetric C-modes 
observed for water rones on Sky lab • 







Fig. 12. Plateau simulation experiments using neutrally buoyant carbon tetrachloride 
plus benzene solution drops in water, (a) and .(c) are static sessile 
drops, (b) and (d) are drops (a) and (c) respectively under vibrational 
conditions. 






Fig. 13. Oscillation frequencies and decay rates of drops such as those in 

Fig. 12. Solid curves theoretically determined by Miller and Striven.™?) 
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